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Abstract 

In the Closest String problem one is given a family S of equal-length strings over some 
fixed alphabet, and the task is to find a string y that minimizes the maximum Hamming dis¬ 
tance between y and a string from S. While polynomial-time approximation schemes (PTASes) 
for this problem are known for a long time [Li et ah; J. ACM’02], no efficient polynomial-time 
approximation scheme (EPTAS) has been proposed so far. In this paper, we prove that the exis¬ 
tence of an EPTAS for Closest String is in fact unlikely, as it would imply that FPT = W[I], 
a highly unexpected collapse in the hierarchy of parameterized complexity classes. Our proof 
also shows that the existence of a PTAS for Closest String with running time /(e) • n°e/e)^ 
for any computable function /, would contradict the Exponential Time Hypothesis. 
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1 Introduction 


Closest String and Closest Substring are two computational problems motivated by ques¬ 
tions in molecular biology connected to identifying functionally similar regions of DNA or RNA 
sequences, as well as by applications in coding theory. In Closest String we are given a family 
S of strings over some fixed alphabet S, each of length L. The task is to find one string y € 
for which max 3 ;g 5 'R(x, y) is minimum possible, where 'H{x,y) is the Hamming distance between x 
and y, that is, the number of positions on which x and y have different letters. We will consider 
both the optimization variant of the problem where the said distance is to be minimized, and the 
decision variant where an upper bound d is given on the input, and the algorithm needs to decide 
whether there exists a string y with max^jg^'H(x, y) < d. Closest Substring is a more general 
problem where the strings from the input family S all have length m > L, and we look for a string 
y E that minimizes max^jg^ min^;/ substring of i 2/)- In other words, we look for y that can 
be fit as close as possible to a substring of length L of each of the input strings from S. 

Both Closest String and Closest Substring, as well as numerous variations on these 
problems, have been studied extensively from the point of view of approximation algorithms. Most 
importantly for us, for both of these problems there are classic results providing polynomial-time 
approximation schemes {PTASes): for every e > 0, it is possible to approximate in polynomial time 
the optimum distance within a multiplicative factor of (l-|-e). The first PTASes for these problems 
were given by Li et al. [9], and they had running time bounded by This was later improved 

by Andoni et al. [T] to and then by Ma and Sun [H] to which constitutes the 

current frontier of knowledge. We refer to the works OEIEIEIIIIKII] for a broad introduction to 
biological applications of Closest String, Closest Substring, and related problems, as well 
as pointers to relevant literature. 

One of the immediate questions stemming from the works of Li et al. [9], Andoni et al. [T], and 
Ma and Sun m, is whether either for Closest String or Closest Substring one can also 
give an efficient polynomial-time approximation scheme (EPTAS), i.e., an approximation scheme 
that for every e > 0 gives a (1 -|- e)-approximation algorithm with running time f(e) ■ for 

some computable function /. In other words, the degree of the polynomial should be independent 
of e, whereas the exponential blow-up (inevitable due to NP-completeness) should happen only 
in the multiplicative constant standing in front of the running time. EPTASes are desirable from 
the point of view of applications, since they provide approximation algorithms that can be useful 
in practice already for relatively small values of e, whereas running times of general PTASes are 
usually prohibitive. 

For the more general Closest Substring problem, this question was answered negatively by 
Marx [12] using the techniques from parameterized complexity. More precisely, Marx considered 
various parameterizations of Closest Substring, and showed that when parameterized by d and 
|5|, the problem remains W[l]-hard even for the binary alphabet. This means that the existence 
of a fixed-parameter algorithm with running time f{d, |<S|) • where n is the total size of the 

input, would imply that FPT = W[l], a highly unexpected collapse in the parameterized complex¬ 
ity. This result shows that, under FPT 7^ W[l], also an FPTAS for Closest Substring can be 
excluded. Indeed, if such an FPTAS existed, then by setting any £< 2 could in time /(d) 
distinguish instances with optimum distance value d from the ones with optimum distance value 
d-|-l, thus solving the decision variant in fixed-parameter tractable (FPT) time. Using more precise 
results about the parameterized hardness of the Clique problem, Marx m showed that, under the 
assumption of Exponential Time Hypothesis {ETH), which states that 3-SAT cannot be solved in 
time 0(2*^"') for some <5 > 0, one even cannot expect PTASes for Closest Substring with running 
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time f{e) ■ for any computable function /. We refer to a survey of Marx [13] for more 

examples of links between parameterized complexity and the design of approximation schemes. 

The methodology used by Marx [T2|, which is the classic connection between parameterized 
complexity and EPTASes that dates back to the work of Bazgan [2| and of Cesati and Trevisan |4| , 
completely breaks down when applied to Closest String. This is because this problem actually 
does admit an FPT algorithm when parameterized by d. An algorithm with running time 
was proposed by Gramm et al. [Tj. Later, Ma and Sun m gave an algorithm with running time 
which is more efficient for constant-size alphabets. Both the algorithms of Gramm 
et al. and of Ma and Sun are known to be essentially optimal under ETH [10], and nowadays they 
constitute textbook examples of advanced branching techniques in parameterized complexity [5]. 
Therefore, in order to settle the question about the existence of an EPTAS for Closest String, 
one should look for a substantial refinement of the currently known techniques. 

An approach for overcoming this issue was recently used by Boucher et al. [3], who attribute 
the original idea to Marx m- Boucher et al. considered a problem called Consensus Patterns, 
which is a variation of Closest Substring where the goal function is the total sum of Hamming 
distances between the center string and best-fitting substrings of the input strings, instead of the 
maximum among these distances. The problem admits a PTAS due to Li et al. [8], and was shown 
by Marx [12] to be fixed-parameter tractable when parameterized by the target distance d. Despite 
the latter result, Boucher et al. [3] managed to prove that the existence of an EPTAS for Consensus 
Patterns would imply that FPT = W[l]. The main idea is to provide a reduction from a W[l]- 
hard problem, such as Clique, where the output target distance d is not bounded by a function of 
the input parameter k (indeed, the existence of such a reduction would prove that FPT = W[l]), but 
the multiplicative gap between the optimum distances yielded for yes- and no-instances is 1 -|- 
for some computable function g. Even though the output parameter is unbounded in terms of k, 
an EPTAS for the problem could be still used to distinguish between ontput instances obtained 
from yes- and no-instances of Clique in FPT time, thus proving that FPT = W[l]. 

Our contribution In this paper we provide a negative answer to the question about the existence 
of an EPTAS for Closest String by proving the following theorem. 

Theorem 1.1. The following assertions hold: 

• Unless FPT = W[l], there is no EPTAS for Closest String over binary alphabet. 

• Unless ETH fails, there is no PTAS for Closest String over binary alphabet with running 

time /(e) • for any computable function f. 

Thus, one should not expect an EPTAS for Closest String, whereas for PTASes there is still 
a room for improvement between the running time of rpU/P) given by Ma and Sun m and the 
lower bound of Theorem ll.il It is worth noting that our /(e) ■ time lower bound for (1 -|-e)- 

approximating CLOSEST String also holds for the more general Closest Substring problem. 
This yields a significantly stronger lower bound than the previous /(e) • lower bound of 

Marx [12]. 

Our proof of Theorem 11.11 follows the methodology proposed Marx m and used by Boucher et 
al. [3] for Consensus Patterns. The following theorem, which is the main technical contribution 
of this work, states formally the properties of onr redaction. 

Theorem 1.2. There is an integer c and an algorithm that, given an instance {G,k) o/C lique, 
works in time 2^ • and outputs an instance {S, L, d) of Closest String over alphabet {0,1} 
with the following properties: 
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• If G contains a clique on k vertices, then there is a string w € {0,1}-^ such that 'H{w,x) < d 
for each x € S. 

• If G does not contain a clique on k vertices, then for each string w G {0,1}-^ there is x ^ S 

such that 'H{w, x) > (1 + • d. 

The statement of Theorem [Olis similar to the core of the hardness proof of Boucher et al. [3]. 
However, our reduction is completely different from the reduction of Boucher et ah, because the 
causes of the computational hardness of of Closest String and Consensus Patterns are quite 
orthogonal to each other. In Consensus Patterns the difficulty lies in picking the right substrings 
of the input strings. Once these substrings are known the center string is easily computed in poly¬ 
nomial time, since we are minimizing the sum of the Hamming distances. In CLOSEST String 
there are no substrings to pick, we just have to find a center string for the given input strings. This 
is a computationally hard task because we are minimizing the maximum of the Hamming distances 
to the center, rather than the sum. 

Theorem o follows immediately by combining Theorem 11.21 with the known parameterized 
hardness results for Clique, gathered in the following theorem, and setting £ = 

Theorem 1.3 (cf. Theorem 13.25 and Corollary 14.23 of [5]). The following assertions hold: 

• Unless FPT = W[l], Clique cannot be solved in time f{k) ■ for any computable func¬ 
tion f. 

• Unless ETH fails. Clique cannot he solved in time f{k) for any computable function f. 

The main idea of the proof of Theorem 11.21 is to encode the n vertices of the given graph G as 
an “almost orthogonal” family 7” of strings from {0,1}^, for some i = O(logn). Strings from 7” are 
used as identifiers of vertices of G, and the fact that they are almost orthogonal means that the 
identifiers of two distinct vertices of G differ on approximately i/2 positions. On the other hand 

1 = O(logra), so the whole space of strings into which V{G) is embedded has size polynomial in n. 
Using these properties, the reduction promised in Theorem ll.2l is designed by a careful construction. 

Notation. By logp we denote the base-2 logarithm of p. For a positive integer n, we denote 
[n] = {1,2,... ,n}. The length of a string x is denoted by |x|. For an alphabet S and two equal- 
length strings x,y over S, the Hamming distance between x and y, denoted 'H{x,y), is the number 
of positions on which x and y have different letters. If S = {0,1} is the binary alphabet, then the 
Hamming weight of a string x over S, denoted H{x), is the number of Is in it. The complement 
of a string x over a binary alphabet, denoted x, is obtained from x by replacing all Os with Is and 
vice versa. Note that if |x| = \y\ = n, then H{x,y) = n — H(x,y) = n — H{x,y) = H{x,y). 

2 Selection gadget 

For the rest of this paper, we fix the following constants: p = 1/100, a = 1/10, /3 = 1/20. Since G is 
divisible by 100, we have that that pi, ai, and jli are all integers. First, we prove that among binary 
strings of logarithmic length one can find a linearly-sized family of “almost orthogonal” strings of 
balanced Hamming weight. The proof is by a simple greedy argument. 

Lemma 2.1. There exist positive integers C and N, where G is divisible by 100, with the following 
property. Let n > N be any integer, and let us denote i = C ■ [logn]. Then there exists a set 
T C {0,1}^ with the following properties: 
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1. \T\=n, 

2. 'H{x) = £/2 for each a: € T, and 


3. (1/2 — p)i < T-i{x,y) < (1/2 + p)i for each distinct x,y €T. 
Moreover, given n, T can be constructed in time polynomial in n. 


Proof. Let H 2 {-) denote the binary entropy, i.e., H 2 {p) = —plogp— (1 — p) log(l — p) for p G (0,1). 
Suppose ^ is some positive integer divisible by 100. Then it is well known that 

k 

^ (A < (1) 

i=0 

for all integers k with 0 < k < £/2; cf. [6l Lemma 16.19]. Let us denote 


A = 


(i/2-p)e 

E 

i=0 


+ 


E 


i=(l/2+p)£ 


Then from ([T]) it follows that 


A<2- 2^^, 


where a = H 2 (l /2 — p) < 1. 

Suppose now that I = C ■ [logn] for some positive integers C and n > 1, where C is divisible 
by 100. Then 


n{i + 1) • ^4 < 

< 
< 


2n-(C'riognl + l)-2‘^'^riognl 

2-{2C + l)- 2^^ • nlogn • 
(4(7 + 2) 


Since fi < 1, we can choose C to 
choose N large enough so that 


be an integer divisible by 100 so that aC + 2 < C. Then, we can 
(4(7 + 2) • 2'^*^ • 


for all integers n > N. Hence, 


nA < 


n 


c 


e + i 


We now verify that this choice of (7, N satisfies the required properties. 

Consider the following greedy procedure performed on {0,1}^. Start with T = 
of {0, lY marked as unused. In consecutive rounds perform the following: 


(2) 

0 and all strings 


1. Pick any x G {0,1}^ with 'H(x) = £/2 that was not yet marked as used, and add x to T- 


2. Mark every y G {0,1}^ with 'H{x,y) < (1/2 — p)i or 'H{x,y) > {1/2 + p)i as used. 

It is clear that at each step of the procedure, the constructed family T satisfies properties (|2]l 
and ([3]). Hence, it suffices to prove that the procedure can be performed for at least n rounds. 

Note that the number of strings marked as used at each round is at most A. On the other hand, 
if P is the set of strings from {0,1}^ that have Hamming weight exactly £/2, then 


V\ > 


1 ( 0 . 1 }' 

t + l 


2Criogn] 

~JTT' - e + i 


From ([2]) we infer that \P\ > nA. This means that the algorithm will be able to find an unmarked 
X gP for at least n rounds, and hence to construct the family T with |T| = n. It is easy to imple¬ 
ment the algorithm in polynomial time using the fact that the size of {0,1}^ is polynomial in n. □ 
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From now on, we adopt the constants C, N given by Lemma 12.11 to the notation. Let us also 
fix n > A^; then let = C • [logn] and T be the set of strings given by Lemma l 2.11 which we shall 
call selection strings. We define the set of forbidden strings T = J'i'T) as follows: 

J" = {y: y € { 0 , 1 }^ and Td{x,y) < (1 — a)i for all x G T}. 

In other words, J- comprises all the strings that are not almost diametrically opposite to some 
string from T. The following lemma asserts the properties of T and J- that we shall need later on. 

Lemma 2.2. Suppose u € {0,1}^. Then the following assertions hold: 

1. If u ^ T, then y) < (1 — a)i for each y G T. 

2. If 'H{x., u) > fli for all x gT, then there exists y G F such that ^{{u, y) > (1 — I3)£. 

Proof. Property ([T]) follows directly from the definition of F, so we proceed to the proof of ([ 2 |). 

Suppose Hiu, x) > j5£ for all x gT. \iu G F, then we could take y = u, so suppose that u ^ F. 
This means that there exists xq G T, for which 'H(xo,Il) > (1 — a)£] equivalently, 'H{'xo,u) < at. 
On the other hand, we have that 'H{xo,u) > /3i, so also TLfxQ,!!) > fli. Construct y from u by 
taking any set of positions X of size /3i on which u and xq have the same letters, and flipping the 
letters on these positions (replacing Os with Is and vice versa). Such a set of positions always exists 
because a + /? < 1. Then we have that Tiixo, y) = 'H{'xo,u) + /3i, which implies that 

a£ = /3£ + (3i < TLixo, y) < {a + (3)£. 

We claim that y G F] suppose otherwise. Since 'H{xo,y) > a£, then also 'H{xo,y) < (1 — a)£. 
As y ^ F, there must exist some xi G T., xq ^ xi, such that ?^(xi,y) > (1 — a)^; equivalently 
'H(xT, y) < a£. Hence, from the triangle inequality we infer that 

7^(xo,xi) = H(xh,xl) < 'H{x^,y) +U{y,xl) < (2a +/3)£. 

This is a contradiction with the assumption that 'H{xo,xi) > (1/2 — p)£, which is implied by 
xo, xi G T. Indeed, we have that 2a + /3 = | < ^ = ll2 — p. 

Hence y G F. By definition we have that 'H{u,y) = ll£, which implies that Ti{u,y) = (1 — I3)£. 
Thus, y satisfies the required properties. □ 

3 Main construction 

In this section we provide the proof of Theorem 11.21 Let (G, k) be the input instance of Clique, 
and let n = |H(G)|. Let C,N be the constants given by Lemma l2.11 We can assume that n > N, 
because otherwise the instance {G,k) can be solved in constant time. Let £ = Cflogn]. We run 
the polynomial-time algorithm given by Lemma l2.II that computes the set T C {0,1}^ of selection 
strings. Let F = F{T) be the set of forbidden strings, as defined in Section [2j Note that F can be 
computed in polynomial time directly from the definition, due to |{ 0 , 1 }^| = 

We now present the construction of the output instance (5, L, d) of Closest String. Set 
L = k£ + 7 ^, where 7 = p + a = and partition the set [L] of positions in strings of length L 

into k + 1 blocks: 

• k blocks Bi for i G [k] of length £ each, where Bi = {{i — 1)£ -|- 1, (i — 1)£ + 2,..., i£}-, 

• special balancing block C of length 7 ^, where C = {k£ + l,k£ + 2,..., L}. 
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For w G {0,1}^ and a contiguous subset of positions X, by w[X] we denote the substring of w 
formed by positions from X. 

Let us first discuss the intuition. The choice the solution string makes on consecutive blocks 
Bi will encode a selection of a /c-tuple of vertices in G. Vertices of G will be mapped one-to-one 
to strings from T. The family of constraint strings S will consist of two subfamilies 5gei and 5adj 
with the following roles: 

• Strings from <Ssei ensure that on each block Bi, the solution picks a substring that is close to 
some element of T- The selection of this element encodes the choice of the ith vertex from 
the fe-tuple. 

• Strings from 5adj verify that vertices of the chosen /c-tuple are pairwise different and adjacent, 
and hence they form a clique. 

A small technical caveat is that for strings from 5sei and from <Sadj, the intended Hamming distance 
from the solution string will be slightly different. The role of the balancing block G is to equalize 
this distance by a simple additional construction. 

We proceed to the formal description. Since |V(G)| = |T|, let t: V{G) —T be an arbitrary 
bijection. 

The family 5sei consists of strings a{i,y,(j), z), for all i G [k], y G X, (f) being a function from 
[A;]\{i} to {0,1}, and z being a binary string of length 7 ^. String a(i, y, cj), z) is constructed as follows: 

• On block Bi put the string y. 

• For each j G [k] \ {i}, on block Bj put a string consisting of i zeroes if 4>{j) = 0, and a string 
consisting of i ones if (^{j) = 1 . 

• On balancing block G put the string z. 

Thus, |<Ssei| = /c • |W| • ■ 2 '*'^ < 2 ^ • Also, 5sei can be constructed in time 2 ^ • directly 

from the definition. 

The family 5adj consists of strings b{i, j, {u,v),'tp), for all i,j G [k] with i < j, {u,v) being an 
ordered pair of vertices of G that are either equal or non-adjacent, and ip being a function from 
[^] \ {b j} to {0) !}• String b{i,j, {u,v),ip) is constructed as follows: 

• On block Bi put the string t(u). 

• On block Bj put the string i{v). 

• On block Bq, for g € [fe] \ {i,j}, put a string consisting of i zeroes if ip{q) = 0, and a string 
consisting of £ ones if ip{q) = 1 . 

• On balancing block C put a string consisting of zeroes. 

Thus, |5adj| < ( 2 ) ■ ■ 2^“^ < 2^ • . Again, 5adj can be constructed in time 2^ • directly 

from the definition. 

Set S = 5gei U<Sadj and d = {k/2 -|- 1/2 + p) G. This concludes the construction. Its correctness 
will be verified in two lemmas that mirror the properties listed in Theorem 11.21 

Lemma 3.1. If G contains a clique on k vertices, then there exists a string w G {0,1}^ such that 
'H{w, x) < d for each x G S. 
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Proof. Let {ci,C 2 ,... ,Ck} be a /c-clique in G. Construct w by putting t(cj) on block Bi, for each 
i £ [k], and zeroes on all the positions of the balancing block C. 

First, take any string a = a{i,y,(j),z) £ iSsei- Since i(cj) £ T and y £ F,hy Lemma [T 2 tfT]l we 
infer that 'H{w[Bi],a[Bi]) = 'H{i{ci),y) < (1 — a)l. For each j £ [k] \ {z}, since = ^/2 due 

to t(cj) £ T, we have that 'H{w[Bj\,a[Bj\) = 'H{i{cj),a[Bj\) = £/2, regardless of the value of 
Finally, obviously T-L{w[C],a[C]) < IC] = 'y£. Hence 

T-Liw, a) < (1 — a)i + {k — l)i/2 + ji = d. 


Second, take any string b = b{i,j, {u,v),'iIj) £ <Sadj- Since q and cj are different and adjacent, 
whereas u and v are either equal or non-adjacent, we have (cj, cj) ^ {u, v). Without loss of general¬ 
ity suppose that c* 7 ^ u] the second case will be symmetric. Then 'H{w[Bi], b[Bi]) = i{u)) < 

(1/2 -|- / 9 )f, due to property ([3]) of Lemma [2Tl Obviously, T-L{w[Bj\.,b[Bj]) < \Bj\ < 1. Finally, for 
every q £ [A:] \{i, j} we have that 'H{i{cq)) = ^/2, and hence 'H{w[Bq], b[Bq]) = T-L{L{cq),b[Bq]) = £12., 
regardless of the value of Strings w and b match on positions of C, so T-L{'w[C],b[C]) = 0. 

Summarizing, 

n{w, b) < (1/2 + p)£ + £ + {k- 2)t/2 = d. 


□ 


Lemma 3.2. If there is a string w £ {0,1}^ such that 'H{w,x) < d + ^£ for each x £ S, then G 
contains a clique on k vertices. 

Proof. We first prove that on each block Bi, w is close to selecting an element of T. 

Claim 1. For each i £ [k\ there exists a unique Xi £ T such that PL{w[Bi],Xi) < fl£. 

Proof. Uniqueness follows directly from property ([3]) of Lemma 12.11 and the triangle inequality, so 
it suffices to prove existence. 

Let u = w[Bi]. For the sake of contradiction, suppose PL{u,x) > j3£ for each x £ T. From 
Lemma [T2l![^ we infer that there exists y £ F such that 'H{u, y) > {l—j3)£. Let us take (/: [A:]\{f} —>■ 
{ 0 , 1 } defined as follows: 4>{j) = 0 if in re the majority of positions of Bj contain a one, and (p{j) = 1 
otherwise. Also, define 2 ; = w[G\. Consider string a = a{i,y,(f,z) £ 5sei. Then, it follows that 

• n{w[Bi],a[Bi]) = n{u, y)> (I- I3)£; 

• 'H{w[Bj], a[Bj]) > £/2 for each j £ [k] \ {i}; 

. n{w[G],a[G]) = niw[G],^JJ\C]) = \G\ = 7 L 

Consequently, 

H{w, a) > (1 — f3)£ + {k — l)£/2 + ^£ = d + (3£. 

This is a contradiction with the assumption that 'H(rc, x) < d + fil for each x £ S. j 

For each i £[k], let Cj = L~^{xi). 

Claim 2. For all i,j £ [k] with i < j, vertices Ci and Cj are different and adjacent. 

Proof. For the sake of contradiction, suppose c* and Cj are either equal or non-adjacent. Define 
V’- [k]\ {i,j} — ^ { 0 , 1 } as follows: ijj{q) = 0 if in rc the majority of positions of Bq contain a one, 
and f>{q) = 1 otherwise. Then, for {ci,Cj) we have constructed string b = b{i,j,{ci,Cj),'tp) £ 5adj- 
Observe now that 
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• 'H{w[Bi],h[Bi]) ='H{w[Bi],Xi) > (1 -/3)£, since 'H{w[Bi],Xi) </3£; 

• Similarly, 'H{w[Bj],b[Bj]) > (1 — I3)i; 

• 'H{w[Bg],b[Bg]) > £/2 for each g € [A:] \ 

• n{w[c],b[c]) > 0 . 

Consequently, 


niw, b) > 2(1 - p)i + {k- 2)£/2 = {k/2 + 1 - 2/3)£ > d + /?£. 

This is a contradiction with the assumption that 'H{w, x) < d + f3i for each x £ S. j 

Claim [2] asserts that, indeed, {ci, C 2 ,..., Ck} is a /c-clique in G. □ 

Lemmas 13. II and I, '1.21 conclude the proof of Theorem ll.21 where c can be taken to be any constant 
larger than < | = 40. 

4 Conclusions 

In this paper we have proved that Closest String does not have an EPTAS under the assumption 
of FPT 7 ^ W[l]. Moreover, under the stronger assumption of the Exponential Time Hypothesis, 
one can also exclude PTASes with running time /(e) • for any computable function /. How¬ 

ever, the fastest currently known approximation scheme for Closest String has running time 
) jYl]. This leaves a significant gap between the known upper and lower bounds. Despite 
efforts, we were unable to close this gap, and hence we leave it as an open problem. 
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